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We construct 3-D solutions of Maxwell's equations that 
describe Gaussian light beams focused by a strong lens. We 
investigate the interaction of such beams with single atoms 
in free space and the interplay between angular and quantum 
properties of the scattered radiation. We compare the exact 
results with those obtained with paraxial light beams and 
from a standard input-output formalism. We put our results 
in the context of quantum information processing with single 
atoms. 



I. MOTIVATION 

The ability to manipulate small quantum systems indi- 
vidually is a necessary requirement for quantum comput- 
ing and quantum communication. For example, in order 
to perform single-qubit operations on a particular ion in 
an ion trap quantum computer Jl|-^] one has to focus a 
laser beam to the position of that ion with a sufficiently 
high spatial resolution [Q . Similarly, quantum communi- 
cation protocols using atoms trapped inside optical cavi- 
ties H|J require single-bit and two-bit operations on spe- 
cific atoms. One question that arises is whether strong 
focusing has an undesired side effect, namely, that the 
scattered light contains information about the state of 
the qubit. The fear would be that the laser intensity 
would have to be turned down so much, that the absence 
of a photon from the laser beam becomes in principle 
detectable. 

Conversely, if an atom in free space would indeed be 
able to modify appreciably the state of a light field, then 
this effect could be used to our advantage: a single atom 
could be used to perform quantum-logic operations on 
single photons in free space R . An atom inside an optical 
cavity strongly coupled to a cavity mode is known to be 
able to perform such tasks j8j, but such an experiment 
would be more straightforward to conduct in free space. 

At first sight the prospects of such an experiment seem 
good: the scattering cross section of a two-level atom is 
(7 = 3A 2 /(27r) for light of wavelength A [||, and thus fo- 
cusing light down to an area A < a should be sufficient 
to induce a strong coupling. However, this picture is too 
simplistic. Light beams are transversely polarized, which 
implies that only part of the light entering the interaction 
region will carry the polarization that the atom is sensi- 
tive to. Using a different picture, since the atom would 
emit a dipole pattern (in a given m — > ml transition), 
it would be most efficiently excited by a field matching 



an "incoming" dipole field, as follows from time-reversal 
symmetry. Since a focused laser beam does not have a 
large overlap with such a dipole pattern, the effective 
absorption cross section is smaller than indicated by a. 

Early experiments p0pl| ] on the detection of single 
atom fluorescence did not reach the strong focusing limit 
of A ~ A 2 . Recently, however, impressive progress has 
been made in experiments on single molecules in con- 
densed matter and efficient detection of the fluorescence 
light has become possible |l2]-|l4| . Experiments on single 
atoms aiming at the strong focusing limit are underway 
as well fT|. 

In a recent paper [|l6| we gave the results of explicit 
calculations on the behavior of single atoms in free space 
irradiated by tightly focused light beams. We were par- 
ticularly interested in the quantum aspects of the scat- 
tered light, and in evaluating how much a single atom 
is able to modify the intensity and phase properties of 
the incident light. Here we present all the details of that 
calculation and discuss possible extensions. These re- 
sults will be compared with similar calculations using 
(i) paraxial Gaussian beams | fl7| and (ii) a well-known 
quantum-optical input-output formalism that was used 
in Refs. |LS|]l9f | to study photon statistics and intensity 
correlations of the field emitted by an atom in free space. 
The latter model is basically a quasi one-dimensional 
model, with one spatial variable describing the propa- 
gation of the light beam, and with one additional pa- 
rameter describing the solid angle subtended by the laser 
beam at the atom's position (coinciding with the focal 
point of the light beam). As we will demonstrate, how- 
ever, neither this model nor, as expected, paraxial beams 
accurately represent the case of a strongly focused light 
beam. 



II. STRONGLY FOCUSED GAUSSIAN BEAMS 

Here we wish to calculate the field that one obtains by 
focusing a monochromatic Gaussian (paraxial) beam by 
an ideal strong lens. We do that by expanding the out- 
going field (i.e. the field after the lens) in a complete set 
of modes. In principle one can use any complete set of 
modes to describe exact solutions of Maxwell's equations. 
In view of the cylindrical symmetry of the problem we are 
interested in, the most convenient is to choose a set that 
takes a simple form in cylindrical coordinates. In par- 
ticular, we will use a set of eigenmodes of 4 commuting 
operators corresponding to the following 4 physical quan- 
tities: energy [with eigenvalue hkc = Hu> per photon [Effl], 



1 



angular momentum in the z direction [mh], momentum 
in the z direction [hk z ], and helicity [shk z /k]. The modes 
are thus characterized by the 4 numbers v = [k,k z ,m, s), 
which, once the field has been quantized, play the role of 
quantum numbers. These modes were constructed in j pjj ] 
to clarify the meaning of orbital angular momentum of 
light p2[ , and thus by construction possess simple prop- 
erties under rotations around the propagation (z) direc- 
tion. The complete orthogonal set of modes F v is defined 
such that the free electric field (i.e. the solution of the 
source-free Maxwell equations) can be expanded in this 
set as 



E = 2Re[^ a v F v exp(— iwt))], 



(1) 



with arbitrary complex amplitudes a v . This requires the 
mode functions to be transverse, i.e. V • F v — 0. The 
summation over v is a short-hand notation for 



5>/d fc /dfc z ££. 



(2) 



The dimensionless mode functions F v in cylindrical co- 
ordinates (p, z, (j>) are defined by PH 

P v (p,z,<f>) = ^-^^G(k,k z ,m+l)e.+ 



47T k 

1 sk + k, 
H £ 

477 k 



G(fc, k z ,m — l)e+ 



47T K 



(3) 



where fc t = y'fe 2 — fef is the transverse part of the wave 
vector, e± = (x±iy)/\/2 are the two circular polarization 
vectors, and 

G(k,k z ,m) = J m (hp) exp(ik z z) exp(im(f)), (4) 

with J m the m-th order Bessel function. The mode func- 
tions satisfy the orthogonality relations 



dVF*(f) ■ F v ,{r) = S(k - k')S(k z - k' z )6 mm ,5 ss , /k, (5) 



where the integration extends over all space. This fol- 
lows directly from the fact that the mode functions are 
eigenfunctions of commuting hcrmitian operators. 

For the remainder of this Section, we will consider only 
monochromatic beams with a fixed value of k = 2tt/X. 
For convenience we take kt as a mode number instead of 
k z and we denote the reduced set of mode numbers by 
p = (kt, rn, s), and introduce the notation 



£=/«*££■ 



(6) 



For fixed k the modes F^ are orthogonal in planes 
z ^constant: 



dSF;(f)-F fi ,(r) = 

z— constant 

S(k t - k' t )S mm ,S ss/ /(271-fct), 



(7) 



which is a useful relation for defining the action on light 
beams of an ideal lens positioned in a plane z=constant. 



A. Focusing with an ideal lens 

The action of the lens is modeled here by assuming that 
the field distribution of the incoming field is multiplied 
by a local phase factor 



p = exp(-ikp 2 /2f), 



(8) 



with / the focal length of the lens p3 24 1. An ideal 
parabolic lens would be represented by a phase factor 
(fip — exp(iky/ p 2 + (/ — p 2 /2/) 2 ). In the paraxial limit 
/ 3> p this factor becomes equivalent to <p. It may be 
that the simple lens factor used here, <p, does not give rise 
to the strongest possible focusing, and that tp p would im- 
prove on this. Moreover, actual lens systems designed to 
focus light down to A ~ A 2 consist of multiple lenses (for 
instance, 12 in ongoing experiments |^q| ), partly to ac- 
commodate for the finite size, finite thickness and other 
imperfections of real lenses as compared to ideal lenses. 
We nevertheless for convenience chose a single lens fac- 
tor ip: it allows for analytical evaluations and the class 
of light beams thus constructed does reach the focusing 
limit of A rj a (for instance, see the plots corresponding 
to / = 100A). This is sufficient for our purposes of show- 
ing that even focusing down to an area of the size of the 
atomic cross section does not quite (by about a factor of 
~ 5) lead to the strong effects one may have expected or 
hoped for. 

If in the plane of the lens, say z — 0, the incoming 
beam is given by 



F in = F Q (p, 0), 
then the output field is given by 



^out(r) 



with 



= 2nk t 



kp 2 



z=0 



dSex P [- l ^)F a -F* u . 



(9) 



(10) 



(11) 



This definition is such that the limit of / — ► oo corre- 
sponds to free-space propagation, as follows from the or- 
thogonality relation (0). Note that the field distribution 
( |l0| ) is an exact solution of Maxwell's equations, irrespec- 
tive of the choice for Fq (in particular, we can take the 
incoming beam to be paraxial). 
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If we approximate the incoming beam by a circularly 
polarized (lowest-order) Gaussian beam with Rayleigh 
range z- m with kz- m 3> 1 byjj] its dimensionless amplitude 



kp 2 ' 

F {p,4>) = exp ( -— ) e 4 



(12) 



then Kfj, is given by 

k m = dminkt kz ~^ S J dpp,Ja(k t p) x 

/ kp 2 \ ( kp 2 
exp (-,_ j exp (-_ 

This integral can be evaluated using 



(13) 



dxa;Jo(/?x) exp(— ax 2 ) = — cxp(— /3 2 /4a), (14) 
o 2a 



and gives the result 

kt k z + sk 



(15) 



with 



^0 



zr - iz 

4 + / 2 ' 
4 + / 



2 ' 



(16) 



The delta function <5 m i expresses the fact that a lens 
cannot absorb angular momentum from a cylindrically 
symmetric light beam upon normal incidence |^,^| : the 
index m of the outgoing beam is 1, because the incoming 
beam has one unit of "spin" angular momentum. 

When the paraxial limit is valid for the outgoing beam, 
i.e., when kzR 3> 1, zr and zq correspond, as we will 
show, to the Rayleigh range and the position of the focal 
plane of the outgoing beam, respectively. But also out- 
side the paraxial limit, the focused light beam is char- 
acterized by the two parameters zr and zq. The largest 
component of the output field ( |Io|) is the e + component, 
which is given by 



zr ~ izq 



dk f 



k, 2k 2 



k\ 



k 2 



k 2 



(17) 



x J (k t p)exp ( --^{zr - izo) ) exp(ik z z), 



with F + = F ■ e* . Note here the similarity between Qlj 
and the expression given in p3| for a class of light beams 



1 Here we assume for simplicity that the focal plane of the 
incoming beam and the plane of the lens coincide. 



generalizing Laguerre-Gaussian (LG) beams [jl7|. The z 
and e_ components of the output field are proportional 
to higher-order Bessel functions, J\(k t p) and Ji{ktp), re- 
spectively, and will therefore vanish on the z axis. As we 
will be interested in the interaction of an atom on axis 
with the focused light beams, we will not consider these 
components here, but they may be important in other cir- 
cumstances. Note as an aside that these two components 
represent beams with 1 and 2 units of orbital angular 
momentum, and and —1 units of spin angular momen- 
tum, respectively, so that the total angular momentum 
of the outgoing beam is indeed Ti per photon. See [^9| for 
a discussion how these different forms of angular momen- 
tum are transferred to the internal and external angular 
momenta of an atom. 

Returning to the e+ component, when the paraxial ap- 
proximation is valid for the outgoing beam, we may take 
out a factor exp(ikz), use fc — k z « k 2 /2k, and extend 
the integration limit in ( |l7j ) to infinity. Defining 



zr + i(z - Z ), 



(18) 



these approximations lead to 

zr - iz Q 



F, 



exp(ikz) / dk t 2 



xj (k t p)exp - 



2k 



(zr — izo) exp(z/cz) 



■ exp 



kp 2 

' 1Z; 



■ (19) 



which, as announced, represents a Gaussian beam with 
Rayleigh range zr and its focal plane located at z = 
zq. We can in fact rewrite the exact result dl7j ) into a 
different form that explicitly displays the corrections to 
the paraxial approximation, 



„ / , s Zr — IZn , 

F + = exp(zfcz) [Ft +F 2 - F 3 ], 



(20) 



with 

Fi = 



■2 f 1 kp 2 



exp 



-kp 1 



2z u 



U 91,2 _ l2 

dh- r^-Joyktp) exp ' 



k 2 



exp ( i(k z -k+ -±)z ) - 1 



k 2 7 



F 3 = 



kt 2fc kf 



J (fct/o)exp 



2k 



2k 



(21) 



Outside the paraxial limit, when kzR is not large, the 
focal plane is no longer at z = z but moves towards the 
lens by several wavelengths, as shown in Figures 1 and 
2, where several examples of intensity profiles of focused 
beams are plotted. Furthermore, unlike in the paraxial 
approximation, the shape of the field is not just deter- 
mined by the value of zr, but depends on zq as well. 
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position 

FIG. 1. (a) Field strength \F+\ of strongly focused Gaus- 
sian beams on the z axis as a function of Z = (z ~ zn)/A. Note 
that the maximum field strength of the incoming field Fo is 1. 
The lens is located at z = and is characterized by / = 100A, 
so that zq « 100A. The incoming Gaussian beam has increas- 
ing values of z in /X = 1 x 10 3 , 3 x 10 3 , 1 x 10 4 . . . 3 x 10 5 , resp., 
for the bottom to top curves. This implies for the outgoing 
beam decreasing values of zr ~ 10A, 10A/3, A, . . . A/30, (b) 
Field strength in the focal plane as a function of the trans- 
verse coordinate p/A. 

The plots of the transverse mode profiles show that 
beyond a certain point the width of the field no longer 
decreases with stronger focusing. One cannot focus down 
a laser field to below a certain limit, roughly about half 
a wavelength, with the ideal lens with lens factor (||), 
no matter how small zr becomes. Moreover, one notes 
the asymmetry of the outgoing beam around the focal 
plane, in contrast to a paraxial beam which is symmetric 
in its focal plane. Since there is no a priori symmetry 
under reflections in the focal plane, this fact should not 
be surprising. 



Let us note here that strongly focused light beams in 
general will display phase singularities (rings in space 
where certain components of E vanish), of the kind 
investigated (both theoretically and experimentally) in 
Refs. 30-33 . For Gaussian illumination of a spherical 
lens, however, no such singularities appear Q]. 

Finally, it may also be interesting to consider tightly 
focused donut beams, i.e., beams of light produced by 
focusing an incoming higher-order LG beam |fr7f . In the 
case of an incoming first-order LG beam, the incoming 
field distribution Fq can be written as (again assuming 
its focal plane coincides with the plane of the lens) 



For the coefficients we find then 



(22) 



k z + sk 



f°° P 2 

/ dp — Ji(k t p) x 



exp 



kp 2 \ ( kp 2 



(23) 



where <5 m 2,o indicates that for the + sign in (22) we get 
S m 2 and for the — sign d m o- These delta functions again 
express conservation of angular momentum: the outgo- 
ing beam possesses the same angular momentum as the 
incoming beam |2(||27]], which has one unit of "spin" an- 
gular momentum and ±1 units of "orbital" angular mo- 
mentum. The integral can be evaluated using 



dxx 2 Ji(fix) exp (—ax 2 



JL 

4a 2 



cxp(-/3 2 /4a), (24) 



and gives the result 



kf k z 



S k e 



■ exp 



2k l 



(25) 
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FIG. 2. As Figure 1 but for / 
Zi n /A = 3 x 10 4 ,1 x 10 5 . 



500A and 



For k~ oc 5 m o there is a nonzero field on axis of a different 
polarization than the incoming field: the z component, 
which would be neglected in the paraxial limit is in fact 
the only nonvanishing component on axis. In this case, 
however, the field on axis will be sensitive both to de- 
viations of the lens from an ideal spherical lens, and to 
deviations of the incoming beam from a pure donut beam, 
even in the paraxial limit. For some explicit examples of 
this sensitivity see, for instance, Refs. |2f^,|3^|. 



III. SCATTERING LIGHT OFF OF A SINGLE 
ATOM IN FREE SPACE 

In this Section we will investigate the response of an 
atom located in the focal region of a strongly focused 
laser beam of the form ( |l7| ) at f = fo- We consider a 



Ja 







J e = 1 transition in the atom, as it is the 



simplest case where all three polarization components of 
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the light in principle play a role. For simplicity we will 
assume the atom to be located on the z axis, so that it 
in fact interacts only with a single (e+) polarization com- 
ponent; as mentioned above, the other two polarization 
components vanish on axis. 

We are mostly interested in calculating the second- 
order correlation function for the light field as a function 
of position and time. For that purpose, the Heisenberg 
picture is the most convenient. In the Heisenberg picture, 
the electric field operator can be written as the sum of a 
"free" part and a "source" part p3], 



E = Ef + E s , 
where the free part is given by 

Ef(r, t) = F v (r)a v exp(— iut) + h. 



(26) 



E 



(+) 
/ 



E 



(-) 



(27) 



Here we separated the field in positive- and negative- 
frequency parts and used the spatial mode functions 
F JJ (r) from Eq. (|§|), with a v the annihilation operator 
for mode v. The source part for the case of a J g = — > 
J e = 1 transition is given by |34| 



E { s +) (r) = J2Mr')^(t - \r'\/c), 



(28) 



where ff = f — f*o , and is the atomic lowering opera- 
tor, and the sum is over three independent polarization 
directions % = ±1,0. Equation (^8|) is valid in the far 
field, with ijji(r) the dipole field 



ipii'r) 



ck _ (dj ■ r)r 



(29) 



Here ujq is the atomic resonance frequency, and di = dv,{ 
is the dipole moment between the ground state \g) and 
the excited state |e,-) in terms of the standard unit circu- 
lar vectors, 



U-i = £- 

u = z, 
U\ = —e 



+ - 



(30) 



and the reduced atomic dipole matrix element d. 

Expressions containing the electric field in time- 
ordered and normal-ordered form (as measured using 
standard photon detectors), such as the intensity and the 
second-order intensity correlation function, can be trans- 
formed into what Ref. Q| denotes as O-ordered form, 
where £?i +) is placed to the left of £?j +) , E^ to the left 

of Ei , and where the source parts are time-ordered. 
For instance, if we assume the initial state of the light 
field to be a coherent state then the normally ordered 
intensity can be written as 



I(t,f) = {E { -\t,F) ■ E i+) (t,f)) 

= • M^ZiU) + |«| 2 |Fout(r)| 2 

»,j=-l,0,l 



i=-l,0,l 



(31) 



where a determines the amplitude of the coherent state, 
such that = <xF out . In Eq. @ we introduced 

the retarded time t r = t — |r"|/c, and a l eg = (cr^~) and 
a ee — { a t a j) are expectation values of the correspond- 
ing atomic operators. The three terms in (|3ll ) corre- 
spond to the intensity Id of the dipole field, II of the 
incoming laser beam, and the interference term. Sim- 
ilarly, the second-order correlation function (where we 
now suppress the dependence of the fields on f) 

G^(t,T,r) = 
{E[-\t)E^\t + T)E^{t + r)E[-\t)) (32) 

i,m— x,y ,z 

consists of 16 terms. For r = 0, 7 of those vanish identi- 
cally, and the remaining ones are 

G^(t,0,r) = |a| 4 |F out | 4 + ^2|«| 2 |F out |V* ■ ^liU) 
+ ]T4Re(a* exp^o*))^ • ^\a\ 2 \F OVLt \ 2 a l eg {t r )) 

i 

+ 2|a| 2 (Fout • ^OCt • $i<%(tr)). (33) 

For the evaluation of the atomic quantities we assume the 
atom reaches a steady state. Define Ci = ad* ■ F out (0)/fi 
and the matrices M. and M.\ :2 by 



y _ 



ac*/r 
r/2 + iA 



2 r/2-iA 
M = (M 1 +M 2 )/(M 1 +M 2 + t) (34) 

with r the decay rate of the excited states , 

and A = ujq — uj a the detuning of the laser field from 
atomic resonance. In the steady state we have then 

1 



99 1 + Tr M ' 



fee = <?ggM 



'eg 



r/2-iA 



exp(— iuiot). 



(36) 
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We are mainly interested in finding the maximum ef- 
fect the atom may have on the outgoing beam. We 
therefore consider the case of weak on-resonance exci- 
tation, i.e., |C| <C r and A = 0. We then calculate 
g ( V(T,r) = G^(T,f)/I 2 (r) at r = 0— there is no de- 
pendence on t in the steady state and for simplicity we 
leave the argument t out — as a function of position in 
the far field. Results are shown in Figures 3-4. The 
distance to the atom is fixed at R = 50A for numerical 
reasons. Note that the angular spectrum does depend on 
the precise value of R. Only in the forward direction do 
the dipole field and a laser field display the same asymp- 
totic behavior. In the forward direction, i.e. on the z 
axis, the laser field turns out to overwhelm the scattered 
field, irrespective of how strongly the light is focused onto 
the atom. This may be compared to a similar result for 
classical scattering from spherical dielectrics with light 
focused down to spot sizes larger than 5 times the size 
of the spheres |^5[ . Hence we find that g^(0,r) w 1 for 
forward scattering, which is in sharp contrast with the 
result from |Q which predicts a large bunching effect 
(i.e., ^> 1) for tight focusing. The figure also shows 
that in a perpendicular direction the dipole field domi- 
nates, so that g( 2 \0,r) = for 4> — > 7r/2 (i.e., there the 
light is almost purely fluorescence light, which is anti- 
bunched (HjH]). g^ reaches a maximum around an- 
gles where the scattered and laser fields are comparable 
in magnitude. The oscillations indicate that <?( 2 )(0,r) is 
very sensitive to the relative phase between the dipole 
and the laser field. In fact, maxima in </ 2 ' appear when 
the free field and the dipole field interfere destructively. 
Indeed, this implies that the total field is smaller than 
the laser field, which implies a photon has just been ab- 
sorbed by the atom. The atom is therefore in its excited 
state, and hence one can expect a fluorescent photon to 
appear soon, thus leading to a strong bunching effect. 

Going from Fig. 3 to 4 corresponds to tighter focusing 
{zr decreases by a factor of 2) and we see that 

1. In the forward direction the ratio of the amounts 
of laser and scattered light decreases (but it's still 
much larger than 1). 

2. The region where g^ reaches its maximum moves 
outward to larger angles <f). 

3. The ratio of the amounts of laser and scattered light 
at <fi = 90° increases by a large amount. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 3. Plot (a) gives the relative intensities of the laser 
field, the dipole field and the total field as a function of the an- 
gle <f>/ir with the z axis (i.e, at position r — [_Rsin <f>, 0, i?cos <j>] 
where we chose R = 50A here and for all further calcula- 
tions. The parameters for the incoming beam and the lens 
are / = 500A and Zi„ = 3 x 10 4 A, so that zr — 8.3A and 
20 = 500A, where we chose A = 852nm, corresponding to the 
D2 transition in Cs, and the atomic dipole moment d adjusted 
so as to give the correct corresponding spontaneous emission 
rate P = 2vr x 2.6MHz for the 6P 3/2 states of Cs. Plot (b) 



gives 



g^ 2 \0, r) as a function of (j>/ir. 




FIG. 4. As Figure £n|, but the parameters for the incoming 
beam and the lens are / = 500A and Zi n = 6 x 10 4 A, so that 
z a = 4.2A and z = 500A. 

We can compare these results with those for a Gaussian 
beam with the same beam parameters. Figures 5 and 6 
show that the 3 conclusions still hold. However, a Gaus- 
sian beam exaggerates the amount of light in the forward 
direction (small <f) at the cost of greatly underestimating 
it for larger angles. This implies that the region where 
g( 2 ) reaches its maximum is moved to smaller angles <f> 
for a paraxial beam. 



G 
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"0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

FIG. 5. As Figure 3, but for a paraxial beam characterized 
by the same beam parameters zr — 8.3A and zq = 500A. 
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FIG. 6. As Figure 4, but for a paraxial beam characterized 
by the same beam parameters zr — 4.2A and Zq = 500A. 

We now focus on forward scattering, and plot in Fig- 
ure 7 the ratio of the intensities of the laser field and 
the dipole field, i.e., K = \Ef\ 2 /\E S \ 2 , in the forward 
direction (0 = 0) as a function of the normalized (di- 
mensionless) beam width w, defined as 



(37) 



The laser field intensity is seen to be much larger than 
the dipole field intensity, by at least a factor of ~ 500. 
For a Gaussian beam, on the other hand, the ratio be- 
comes arbitrarily small for small w. This has immediate 
consequences for the value of g^(0, r) (see Figure 8). 



FIG. 7. The relative intensity of the laser beam to the 
dipole field in the forward direction, K, as a function of 
the beam width parameter w for the case / = 500A. Since 
zr < f/2, the beam width w satisfies w < 6.3. The dashed 
curve corresponds to the exact solution, the solid curve to a 
Gaussian beam. 

For a Gaussian beam, the intensity in the focal region 
is not bounded. In fact, for decreasing values of w, more 
and more energy is concentrated in the focal region, so 
much so that the dipole field will eventually dominate the 
field in the forward direction. In that case, the forward 
direction will display anti-bunching (for w smaller than 
approximately 0.07). Before that, however, <?( 2 ' reaches 
a large maximum at about w « 0.2, namely when the 
dipole and laser fields are comparable in magnitude. For 
larger values of w, the scattered field will be negligible, 
and i?*- 2 -* — * 1, but only after reaching a minimum close 
to zero around w = 0.3. The latter characteristics were 
also found in Ref. |18| . For the exact solutions, however, 
none of these effects is present, and the laser field always 
dominates the dipole field, so that g( 2 ' w 1 for all values 
of the beam width. 
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FIG. 8. g*- 2 ' (0, r)) in the forward direction as a function 
of the beam width parameter w. The solid curve corresponds 
to a Gaussian beam, the dashed curve to the exact solution 
&■ 

Let us finally quantify the effects of focused light 
on atoms in a different way by considering the follow- 
ing. If the atomic dipole is d, then the relevant quan- 
tity determining the excitation probability of an atom is 
\d ■ E^(ro)\ 2 evaluated at the atom's position Fq, while 
the total incoming energy flux is given by JdS\E^\ 2 . 
In contrast to the naive expectation R ~ cr/A, the ac- 
tual ratio R s that determines the fraction of the energy 
incident on the atom which will be scattered is given by 



R s 



3A 2 |d-jg(-)(f )|- 
2tt / dS\E(-)\ 2 



(38) 



This ratio is plotted in Fig. 9 as a function of the width 
w = wr/X for several values of the focal parameter /. 
For smaller / the best achievable ratio increases, as ex- 
pected, but for realistic lens parameters, the optimum 
R s is about 20%. Even for small values of /, the maxi- 
mum scattering ratio does not go beyond 1/2. This can 
be understood by noting that the optimum shape of the 
illuminating field would be a dipole held. Here with light 
coming only from one direction, one may expect R s to 
be at most 1/2. Obviously, with one mirror behind the 
atom, one can improve the scattering ratio R s by a factor 
of 2. And of course, by building an optical cavity around 
the atom the atom-light interaction can be enhanced by 
many orders of magnitude, but that's a different story 

M. 



CC 6.25 




FIG. 9. The scattering ratio R s as a function of the nor- 
malized beam width w = wr/A. The different curves corre- 
spond to //A = 2.5, 5, 10, 25, 50, 100, 250, 500, 1000 resp., for 
top to bottom curves. 



IV. DISCUSSION AND CONCLUSIONS 

We constructed propagating wave solutions of 
Maxwell's equations describing tightly focused laser 
beams. The method we used consisted of expanding the 
outgoing beam in a complete set of solutions and match- 
ing it at the plane of a lens to a given incoming beam. 
The lens was assumed ideal (infinitely thin) and the in- 
coming beam was chosen to be Gaussian. 

We then investigated quantum-statistical properties of 
the light emitted by an atom in free space, when it is 
illuminated by such a beam. Light detected in the for- 
ward direction does not display any bunching, nor anti- 
bunching effects: the field is dominated by the laser 
light, and the normalized second-order intensity corre- 
lation function is practically unity. This may not be 
surprising but is in contrast to results obtained by us- 
ing Gaussian beams and by a standard quantum-optical 
input-output model. Gaussian beams are no longer valid 
approximate solutions under strong focusing conditions, 
and in particular exaggerate the focal intensity by a large 
amount. On the other hand, the input-output formalism 
implicitly assumes that the scattered field propagates in 
the same manner as the incident light beam: in free space 
this would correspond to illumination with a laser field 
whose profile mimics the dipole pattern. Inside a cavity, 
however, the model is expected to apply, as the situa- 
tion there is to a good approximation one-dimensional. 
Indeed, the equations ultimately assume the same form 
as those for an atom coupled to a cavity mode in the 
bad-cavity limit |38f| . 

Although the model of a single ideal lens with a sim- 
ple lens factor of Eq. (||) may not lead to the strongest 
possible focusing |H| , the amount of focusing reached is 
sufficiently strong (focusing areas A less than or equal to 
the absorption cross section a — 3A 2 /(27r) of an atom) to 
conclude that the interaction a focused light beam with 
an atom is not as strong as might be expected on the 
basis of the ratio a/ A (see Fig. 9). Its consequence for 
quantum information processing may be phrased as: in 
free space it is easier and more efficient to use light to 
process quantum information carried by an atom, than 
to use an atom to process quantum information carried 
by photons. 
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